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Abstract 

We derive a semiclassical transport equation for fermions propagating in the presence of a 
CP-violating planar bubble wall at a first order electroweak phase transition. Starting from 
the Kadanoff-Baym (KB) equation for the two-point (Wightman) function we perform an 
expansion in gradients, or equivalently in the Planck constant h. We show that to first order 
in h the KB equations have a spectral solution, which allows for an on-shell description of 
the plasma excitations. The CP-violating force acting on these excitations is found to be 
enhanced by a boost factor in comparison with the 1+1- dimensional case studied in a former 
paper. We find that an identical semiclassical force can be obtained by the WKB method. 
Applications to the MSSM are also mentioned. 
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1 Introduction 



Modeling electroweak baryogenesis (EWBG) |J requires a study of the generation and trans- 
port of CP-violating flows arising from interactions of fermions with the expanding phase 
transition fronts. Indeed, the most prominent problem in the EWBG scenario over the 
past few years has been to find a systematic derivation of the appropriate transport equa- 
tions ||, including the CP- violating sources. This paper is the second in a series dedicated 
to reach this goal. In the first paper (KPSW-I) we have derived a semiclassical trans- 
port equation starting with the Kadanoff-Baym (KB) equations |5| for the Wightman 
out-of-equilibrium two-point function. For simplicity we have restricted ourselves to 1 + 1 
dimensions there. An important tool for the derivation was an expansion in gradients in the 
center-of-mass coordinate. This is a good approximation when the typical de Broglie wave 
length £dB ~ 1/T of a particle in the plasma is small in comparison to the wall thickness, 
^dB <S ^waii- This is the case for example in many supersymmetric models, which typically 
have £ wa n ~ (10 — 20)/T ||. Here we extend the derivation to the more realistic case of 
fermions interacting with a planar phase transition front in 3+1 dimensions. As in || we 
only consider collisionless fermions and we do not account for the gauge degree of freedom. 
Some aspects of collision terms were discussed in and the gauge fields were studied in the 
1+1 dimensional case in ||. 

A crucial step in our derivation is to show that, to second order in gradients or equivalently 
to first order in h, the KB equation for the Wightman function G < written in the mixed 
(Wigner) coordinates admits a spectral solution. This is not a generic feature, and does not 
persist beyond the order ft. Nevertheless, to first order in K one can write an effective kinetic 
equation for a suitably chosen distribution function for the corresponding on-shell excitations. 
This kinetic equation has the familiar Liouville form, but with a semiclassical force term that 
includes a CP- violating part of order h which sources baryogenesis. This CP-violating force 
is in agreement with the earlier results || derived by use of WKB-methods JTD], [TTJ, |T2| . 

It should be stressed that the existence of the spectral solution to order H is a new result. 
Indeed, while fermion dynamics in the presence of classical background fields have been 
studied extensively within the Schwinger-Keldysh formalism [0, |14], |15|, [16|, 17fl , apparently 
no attempts have been made to find explicit solutions beyond the classical approximation. 
Neither the fact that the spectral decomposition ansatz fails beyond the order h seems to 
be known. Our result is also crucial for EWBG-calculations since no CP-violating terms 
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appear in the classical mean field approximation. Moreover, as the failure of the spectral 
decomposition ansatz beyond order h means that higher order effects cannot be described in 
the semiclassical limit, we are very lucky that in the most interesting cases for the EWBG the 
walls are thick enough so that the gradient expansion converges rapidly and our semiclassical 
theory can be expected to be a very good approximation. 

This paper is organized as follows. In section 2 we derive the Kadanoff-Baym equation 
for a collisionless fermionic system and introduce the Lorentz-boost that relates the wall 
rest frame and the local frame with no parallel momentum, where the KB equations can be 
reduced to the 1+1-dimensional form. We then construct an explicit connection between the 
two frames and make use of our old 1+1 dimensional results [[| to derive the semiclassical 
transport equations appropriate for the EWBG-problem in 3+1 dimensions. In section 3 we 
consider the transformation properties of fermionic currents and relate the CP-violation in 
the fluid equations to the non- conservation of the axial current divergence. In section 4 we 
derive the semiclassical force in the WKB-approach and in section 5 we apply our results to 
chargino transport in the MSSM. Finally, section 6 contains our conclusions. 

2 Derivation of kinetic equations 

As in KPSW-I [|J, we will consider a fermionic field described by the Lagrangian 



arises from an interaction with a spatially varying scalar field condensate. This situation 
is realised at the first order electroweak phase transition (EWPT), which proceeds via nu- 
cleation and growth of the broken phase bubbles of a nonzero Higgs condensate. EWPT 
bubbles quickly grow very large in comparison with the thickness of the phase transition 
front, so one can to a good approximation assume a planar symmetry on the scale of the 
bubble walls. As a result we can take m to be a function of one spatial coordinate. More- 
over, bubble walls are often found to be very wide || in comparison with the typical de 
Broglie wavelength of a particle in the plasma, £ wa n <C ~ l/T. This suggests that one 
should be able to treat the kinetics of plasma excitations near the phase transition front in 
an expansion in gradients. 





where the mass 




(2) 
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Here we consider only a collisionless plasma. Interactions are specified by the Lagrangian 
Ci n t, and in addition to reproducing the usual collision terms in kinetic theory, they give rise 
to a new CP-violating source f|. We will study these effects in detail elsewhere [[L8|]. For 
simplicity we also restrict ourselves to a case without gauge degrees of freedom. Transport 
equations with a self-consistent electric field were considered in 1+1- dimensions in Ref. ||. 

As in [[|, the basic object of our study is the dynamical two point (Wightman) function 

Gap{u,v) = i(ipp{v)tp a {u)) , (3) 

where the expectation value is taken over the initial state. In the collisionless case the 
equation of motion for G < can be obtained by use of the Dirac equation that follows from 
the Lagrangian ([[]): 

(i $ u -m R (u) + zm/(tt)7 5 )V , («) = 0. (4) 

In order to establish the gradient expansion, we first separate the average center-of-mass 
coordinate x = (u + v )/2 from the relative coordinate r = u — v by performing the Wigner 
transform (Fourier transform with respect to r) to the mixed representation: 

G < (x,k) = J d*re ik - r G < (x + r/2,x-r/2). (5) 

For further reference we note that G < (x, k) satisfies the hermiticity property 



t 



l ^G < \ = i 1 {i G < (6) 

implied by the definition Using (Q) and transforming to the mixed representation, one 
finds the equation of motion 

-p-m R (x- ~d k ) -i^m^x- -d k )jiG < = 0, (7) 



where 



m Ri i(x - l -d k ) = m RJ (x)e ^ 9 ' 9k 



Formally equation ([?]) is completely general and valid up to any order in gradients and for 
arbitrary space and time dependence of the mass term (§). This expansion in gradients of 
the center-of-mass coordinate x can be also viewed as an expansion in powers of the Planck 
constant. We have set h — > 1, but a dimensionful h can at any stage be easily restored by 
the simple replacements d x — > hd x and G < — > G < . 
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In the EWBG problem the system is both stationary and has planar symmetry along 
spacelike hypersurfaces of constant z in the wall rest frame. This means that in the plasma 
rest frame the mass has the form 

™>r,i = mn,i(lw(z - v w t)), (9) 

where v w is the velocity of the phase transition front and 7^ the corresponding boost factor. 
It is easiest to perform the analysis in the wall rest frame where rriRj = mnj(z). Moreover, 
because equation is Lorentz-covariant, one can immediately write the equations of motion 
in the wall frame. We also assume a stationary solution for G < , so that its functional 
dependence can be parametrized in the wall frame as G < (k fJi ;z). Q Keeping the gradient 
corrections still to arbitrary order the equation of motion (0) then becomes 

(7% - 7 3 (A; 2 - % -d z ) - 7 • %\\ - m R (z - ~d kz ) - i^m^z - ~£> fe j) iG< = 0. (10) 

Despite of these simplifications, equation fllCf ) still appears quite formidable to solve. In- 
deed, the Wightman function G K contains 16 independent real functions and equation ( |10D 
consists of 16 complex equations for these components^. However, when fen — > 0, i.e. for 
particles that move orthogonally to the wall, equation (|T0|) reduces to the much simpler 1+1 
dimensional one studied in || . We can make use of those results in the general case studied 
here if we can find a transform that removes the 7 • k\\ -term from equation flTU|). In the static 
case under consideration this is done by the boost A to the frame in which kii =0. Obviously 
the four momentum in the boosted frame becomes 

fe|| = 0, k z = k z , fc = sign[k ](A; 2 -fc | f)i (11) 

We define the boost velocity v» = k\\/ko and the gamma-factor 711 = fco/^o f° r later reference. 
Finding the spinor representation of the boost A requires a little more work, but it can be 
deduced from the fact that the desired operator L(A) must commute with 7 s and 7 s and it 
must effect the transform 

L(A)^L- 1 (A) = j%-ry 3 k z = % (12) 

The commutation requirements immediately suggest the form L(A) = a + &7°7 • and 
imposing the condition ( |l"2"D then defines the coefficients a and b, leading to the operator 

L(A) = fc ° + ^- 7 °^. (13) 
'2k (k + h) 



1 To what extent the assumption of stationarity can be relaxed is discussed at the end of this section. 
2 To get a flavour of how the complete set of equations looks like, we refer to Refs. |h| [lTy . 
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The inverse boost operator is simply L l (A(ko,k)) = L(A(&o, — k)). After performing the 
boost, we finally arrive at the following effectively 1+1-dimensional equation 

(fh - l\K - l -d z ) - m R (z - ~d kz ) - vfmfe - ~d kz )) iG< = 0, (14) 

where the boosted Wigner function G < is related to the original one by 

G< = L(A)G < L(A)- 1 . (15) 

The further analysis of equation (|TJ]) is identical to our earlier treatment || . However, 
for completeness we review the main physical aspects of the derivation here. The most 
important property of the 1+1-dimensional frame is that the differential operator in (|i~4"D is 
entirely spanned by the 1+1-dimensional Clifford algebra, and that it commutes with the 
operator 

S z = 7 °7 3 7 5 , (16) 

which measures spin s in ^-direction in which the wall moves. Therefore s is a good quantum 
number in the boosted frame, and as a result we can restrict ourselves to finding solutions 
for G < that satisfy S z Gf = GfS z = sGf. This condition immediately leads to the following 
decomposition for G < , which is block-diagonal: 

- l ^G< = \{l + sa i )®p a ~g s a . (17) 

Here we made use of the chiral representation for the gamma matrices: 7 — > p 1 , — 27°7 5 — > 
p 2 , — 7 5 — > p 3 and 5* z = 7°7 3 7 5 — > a 3 , where p a = (l,p l ), and a 3 and p % are the Pauli 
matrices. The normalization was chosen such that the component §q corresponds to the 
phase space density of states of spin s in the boosted frame. All of the results derived in || 
apply for the functions g s a of course. From these we only need the constraint equations (here 
written to first order in gradients) 

^o$S - sk z9 S 3 ~ m R9i ~ m i92 = (18) 

^3 - sk z9o + 2 m 'A<?2 - ^ m 'idkJl = (19) 

s 1 
ko9t + ^d z g s 2 - m R g s + -m'jd^gl = (20) 

s 1 
ko9 s 2 - 2 dz 9t ~ m i9o ~ ^rn' R d k J s 3 = 0. (21) 
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Equations ( |19l - |2Tl) allow us to express the functions gf in terms of g$ alone, while the re- 
maining equation fll8f) contains the spectral condition for g$. 



We now make an explicit connection between the solutions in the boosted and unboosted 
frame. To this end it is important to note that in the stationary case under study the 
spin in z-direction remains a good quantum number also in the unboosted frame. This is 
mathematically expressed by the fact that the spin operator in the original frame 



L- 1 (A)5 2 L(A) = 7 ||(&-i(i7j| xa),) 



(22) 



commutes with the differential operator in equation (|T0|). As a consequence, even in the 
original frame, the problem splits into two non-mixing sectors labeled by the spin s. By 
construction the Wightman function 



Gf = L(A)^GfL(A) 



(23) 



in the original frame commutes with the spin operator S z . Note that the spin S and a = 7°7 
transform under boosts as the magnetic and electric components of an antisymmetric tensor, 
respectively. This can be seen from the relations a % = —icr 0t and Si = (l/2)ejjfcO" jfc , where 
u^u _ (2/2) [7^, 7"] is the totally antisymmetric tensor. 

Now, the most general decomposition for the Wightman function Gf in the unboosted 



frame can be written in the chiral representation as |T3|, [L4|, [16 

1 „ 



ry°G 



-a 



P b 9 S ab, 



(24) 



where p a = (l,p*) and a a = (1,0"') 0. The connection between the component functions 
in (|24|) and fllTD is now founcl^ by direct use of the relation (pl|) and the explicit form of the 
boost operator (O): 



9ab 



( l\\9o 

vxi\\g s z 
v 2i\\gi 

\ sg s o 



91 



92 Ws \ 

sv 2 i\\gj -svaTll^f v H\\9o 

svi-f\\g s 2 sv^pl v 2 7j|#o 

si\\9l S7||^| sgl j 



(25) 



3 It is easy to show that p a <8> <J b span the full 3+1 dimensional Clifford algebra. 

4 An alternative way to derive this connection would be to identify the transformation properties of g s ab 



by considering the fcrmionic currents, see section (H). 



6 



Given the transformation ( p5f) and the constraint equations (18|-2lf) , we are now ready to 
deal with the equations of motion fllOD . Out of the 16 equations we are primarily interested 
in the 00-component, defining the kinetic properties of (Jq , which corresponds to the particle 
number density in the phase space. Other components can always be determined from (7q 
by the use of constraint equations, and their physical significance will be discussed in the 
next section. The relevant kinetic equation has the form 

9 *9a3 - m 'A#oi - m 'i d k z 9o2 = ( 26 ) 
and the constraint equation reads 

Moo - k i9 s i3 - m R9m ~ mi9o2 = °- ( 27 ) 



In ( p6|) we have kept terms to the second and in (p7]) to the first order in gradients, which 
suffices for the derivation of the kinetic and constraint equations accurate to order h. Using 
(p5|) and the constraints (|18| - ^I1 ) it is now a simple matter to show that we get the following 
kinetic equation 

hd z g s 00 - {hmfd kz - -t(H 2 0')'^) = 0, (28) 
and the algebraic constraint equation 

(k 2 -\ m \* + ^(\ m \ 2 6'))g s 00 = 0, (29) 

which both contain nontrivial CP-violating terms oc 9', where the angle 9 corresponds to the 
complex phase of the mass term m = \m\(z)e ie ^ . 

The most important outcome of these results follows from the observation that the equa- 
tion (^J) has a spectral solution^ 

#oo = ^2-z—n a 6(ko^u a ±), (30) 



where to s ± denotes the dispersion relation 



1 777.1 2 $' 

uj s ± = looTs——- — , uj = \/k 2 + \m\ 2 (31) 

2w (^o - ) 1 



5 The normalization of G < in (|l7|), combined with the spectral sum rule ij° J dk( j (G > — G < ) = 2nl, fixes 
the normalization in (^fj|) unambiguously. It differs by a factor of 4 from our definition in |fj] , which was thus 
not entirely consistent. 
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and Z s ± = 1 =F s\m\ 2 9' /2(u!q — k^) 3 ^ 2 . Because of the delta function in (|30| ) the functions 
n s (ko, k z , z) are projected on-shell and become the distribution functions f s+ and f s - for 
particles and antiparticles with spin s, respectively, defined by 

f s+ = n s (uj s+ ,k z ,z) 

/ s _ = 1 - n s (-u s _, -k z ,z). (32) 

This on-shell projection proves the implicit assumption underlying the semiclassical WKB- 
methods, that the plasma can be described as a collection of single-particle excitations with 
a nontrivial space-dependent dispersion relation. 

Integrating (^8|) over the positive and negative frequencies and taking account of ( p0[ ) 
and Q3"2]) we now get the following on-shell kinetic equations: 



v s ±d z f s ± + F s ±d k J s± = 0, (33) 

where the quasiparticle group velocity v s ± = k z /u s ± is expressed in terms of the kinetic 
momentum k z and the quasiparticle energy u s ± fl3T|) , and the semiclassical force 



m\ 2 ' ± s{\m\H')' ^ 



S± 2u s± 2u (co 2 -k 2 )i 

When compared with the 1+1 dimensional case the sole, but significant, difference in the 
force is that the CP- violating ^'-term is enhanced by the boost-factor 7y = uo/^u 2 — k 2 )^ . 

So far we have considered the Boltzmann equation ([33]) for a static situation (in the wall 
frame) with planar symmetry It is interesting to see to what extent these restrictions can be 
relaxed. The essential aspect of our treatment above was the requirement of spin conservation 
in ^-direction, guaranteeing the consistency of the spin conserving decomposition G < — > Gf 
of the Wigner function. This requirement can be formally expressed by the condition 

[S z ,V]Gf = 0, (35) 

where T> is the derivative operator in ([?]), and the commutator evaluates to 

[S z , V] = -7||7°i7j| x d S[l + x (d^ + v\\ d t ). (36) 

The condition (|35| ) is clearly satisfied when the parameter dependencies of the Wigner func- 
tion are of the form f] 

Gf = Gffa z, t - v {] ■ x || ). (37) 

6 We have checked that, when (|37]) is inserted in VGf = 0, one obtains a fully consistent set of 16 
equations for the components g s a (a = 0..3). 
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For static, planar symmetric solutions studied so far the commutation condition is satisfied 
trivially, because for them dtG < = d x ^G < = 0. We note that in the 1 + 1-frame (in which 
fcii = 0) our result (j37|) is particularly simple and intuitive. In this frame di is spin conserving, 
while dg„ violates spin, and thus the only off-diagonal contributions in the equation T>Gf = 
are the dg.. -derivatives, leading to the requirement dg..g% = 0. This is completely consistent 
with fl37D , which in the 1+1-frame reduces to Gf = Gfik^; z,t). 

Following the procedure outlined earlier in this section one can show that there are no 
new corrections to the constraint equation (|29|), implying that the transients ( |37D respect 
the simple algebraic on-shell condition. The kinetic equation (|33| ) now generalises to 

d t fs± + v s ± ■ dgf s ± + F s ±d k J s ± = 0, (38) 

where v s ± = k/u s ±, and the distribution function satisfies f s ± = f s ±(k^; z, t — v\\ -x\\), where 
now v» = k\\/u a ±. We shall assume this to be satisfied when discussing the current continuity 
equations in the next section. 

3 Fermionic currents and CP-violating sources 

In order to shed light on the physical significance of the various components of the matrix 
Gf, it is instructive to consider the fermionic currents 



^(x)r^(x)) = -J ^Tr [TiG<(x, k)) (39) 

T= (l,7 5 ,7 M ,7V,^ = ^[7 M ,7l)- (40) 

By direct computation one sees that and g^ 2 correspond to a scalar and pseudoscalar 
density, respectively, while (<7q , g? 3 ) forms a 4- vector, and (gg 3 , g^) a 4-pseudovector. Finally, 
g^ and gf 2 are components of an antisymmetric tensor, and so they transform like the 
magnetic and electric field, respectively. These transformation properties could actually 
have been used as an alternative and elegant way to infer the relation between the boosted 
and nonboosted frames in the previous section. 

Let us now look at the vector and pseudovector currents in a little more detail. By 
writing g s ab in terms of an d dropping the total derivative terms we find 

d 4 k 



35s = I (a^iMl' ^ sv ^\\9 S oo- (41) 
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These expressions are clearly the expected generalizations with appropriate boost factors of 
the 1+1 dimensional currents derived in ||. 

It is straightforward to show from the equations of motion, that the vector current diver- 
gence is conserved, and the axial current has the usual form (since we do not treat the gauge 
fields here, the axial current is taken to be nonanomalous; for an account of the anomaly 
see i): 

d,f s = (42) 
d^ s = 2im R $ sl 5 ij s )-2mS s il) s ). (43) 

If we define the fluid density[] n s ± and fluid velocity moments as[| 

«-> - — I T^f^oo, (44) 



n s+ J+ (27r) 4 \k 

then Eqs. ( f42"l - fl3"|) can be shown to correspond to the two lowest velocity moments of the 



Boltzmann equation (38 



d t n s ± + d z (n s ±(v s ±)) = 0, (45) 
d t (n s± (v s± )) + d z (n s ±(v 2 s± )) = S a ±. (46) 

The quantity S 3 ± is related to the baryogenesis source and is given by the average over the 
semiclassical force (0) divided by k : 

S s ± = ~\m\ 2 % s± ±^s(\m\ 2 eyi 3s± , (47) 

and 

It should be stressed that while equations (^5|-|46|) are illustrative in relating the baryo- 
genesis source in moment equations to the non-conservation of the axial current, they are 



7 n s ± should not be confused with the phase space density n s in (p0|). 

8 Here we defined integration over the positive frequencies. To obtain equations for antiparticles one 
integrates over negative frequencies taking into account the appropriate thermal equilibrium limit in (pj 
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somewhat formal, and not directly suitable for actual calculations. For practical baryoge- 
nesis calculations it is convenient to treat the Boltzmann equation fl38f ) in linear response 
approximation with respect to deviations from thermal equilibrium 

8fs± = fs±~ fos±, fos± = e f3 lw (uj a± +v w k z ) + ( 49 ) 

where Jq s ± is the equilibrium distribution function in the wall frame, f3 = 1/T, v w is the 
wall velocity and j w is the corresponding boost factor. Working to linear order in v w we can 
rewrite (|38|) as 



(d t + v s ± ■ d s + F s± d kz ) 5f s ± + v w F s ±{d UJ fu)u>=u B ± = 0, (50) 

with f u = l/(e^ w + 1). Equation (|50D implies the following kinetic equation for the CP- 
violating density 5f s = 5f s+ - 



k 



\m 



2' 



dt H dg- — — d kz I 5 f B = -v w ( 5F s {dufu)u Q + F Su s 



'd u f u 



^0 



(51) 



where we kept only the leading second order terms in gradients, and 

12' 



F 
8uj s 



\m\ 



2u 
\m\ 2 6' 

o; o)o 



6F S = F s+ -F s ^ = s^T^, (52) 

where ujq = {uJq — k^) 1 ^ 2 . Equation (^1|) displays two very important features not evident in 
equations (f45| - [46| ): the CP- violating perturbation Sf s is only sourced by terms second order 
in derivatives (first order in ft) and the source is suppressed by the wall velocity. 

To get equations suitable for baryogenesis calculations it would be necessary to include 
collision terms in (0), which are responsible for all thermalisation and transport processes. 
While treating collisions is beyond the scope of this paper, let us nevertheless mention that an 
approximate solution for Eq. flSTp with collisions is typically found by taking moments of the 
kinetic equation, leading to a set of fluid equations for the moments. In order to get a good 
approximation while using only a few lowest moments, it is important to parametrize 8f s such 
that the parametrization models the actual solution in the best possible way. Particularly 
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suitable and often used is the fluid ansatz parametrised by the chemical potential and a 
perturbation that is an odd function of the momentum. The chemical potential then captures 
an excess of particles formed at the phase boundary front, while the latter parametrizes the 
anisotropic bulk motion of the fluid. For a more detailed account of the fluid equations we 
refer the reader to Refs. [0, 0, g, . 



We finally note that the continuity equations for the vector and axial vector current (45 - 



46) are equivalent to the two lowest order moments of the Boltzmann equation (|5T|) 3, and 



hence also to the fluid equations. 

4 WKB method 

For the sake of comparison we shortly review the derivation of the semiclassical transport 
equations via the WKB- method |L0j, [□], |[], in which the dispersion relations of the single 
particle excitations are derived directly from the Dirac equation (|4]) in the gradient expansion. 
In the wall frame the mass (|2]) is only a function of z, m(z) = \m{z)\e l6 ^ z \ and therefore the 
energy and momentum parallel to the wall are conserved quantities. We can then seek the 
solutions of the form ip{x) = exp(—ik t + ik\\ ■ x\\)tp ko j:^(z), which gives 

(7% - 7 • k\\ - i^d z -m R - Z7 5 m 7 ) (z) = 0. (53) 

As in the previous chapter, we can use the Lorentz boost A (|H|) to get rid of the 7 • k\\ term 



in (|53D . The remaining part of the Dirac operator then commutes with the spin operator S z 



(|j~6l), indicating that the spin s in ^-direction is a good quantum number, and one can write 



the boosted wave function ip kQ = L{h)i}}^,. as a direct product of 2 x 2 spinors, 

h=n ®xs, < = ( 54 ) 

where Xs is the eigenspinor of the spin operator (|l6|) in the 2x2 spinor space, cr 3 x s = sx s , 
and k = sign[fc ](fco — k^) 1 ^ 2 as usual. With this decomposition Eq. (^) simplifies to the 
following two equations for the left and right chiral states in the boosted frame: 

(k — isd z )L s = mR s 

(k + isd z )R s = m*L s . (55) 

Our analysis so far differs from the ones in the literature in that we have discussed more 
explicitly the boost and its relation to the spin-decomposition. The following steps in the 
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analysis are quite standard, and we refer the reader for example to f| U for more details. 
Solving R s in terms of L s in (^) and then introducing the WKB-parametrization 



u,e 



one can write down a set of two coupled second order equations for u s and k s . Solving these 
equations to second order in gradients gives the dispersion relation 

(gfco + p )fl' . 

k s =Po + — , (57) 



where po = sign(A; s ) y — \m\ 2 and s^= 1(— 1) for particles (antiparticles). Remembering 
/cq = /cq — k 2 and inverting fl57| ) one can write the energy in the unboosted frame as 




^ \/ (K - s dp — ) 2 + k\ + \m\ 2 - s^— . (5i 



2 

From this expression for the energy we can finally compute the physical momentum, which 
is defined in terms of the group velocity v g = ko(dk s k Q ) z : 

s\ui\^8' 

k z = k v g =p + - (59) 

7 2p k 

and the corresponding physical force F z = k z = u) s ±dv s ±/ dt: 

\m\ 2 ' s(|m| 2 ^)' , x 

2k <P 2k k 1 ; 

These results are in complete agreement with those found in section 2, completing the explicit 
proof of correctness of the semiclassical limit in the present paper. 

Let us finally point out that the kinetic equations could be formulated in terms of the 



canonical momentum (57), as it was done originally in [[11], 0. The transformation that 



relates the kinetic and canonical momenta is noncanonical: 

kz = hs { lT 2(fcg + H 2 )^ ) • (61) 

The advantage of the canonical formulation is that it allows for a hamiltonian formulation of 
kinetic equations, while the disadvantage is that the canonical momentum and energy do not 
transform as a 4- vector under Lorentz boosts. This makes the formulation of fluid equations 



cumbersome, and so far it has not been done with due care in the literature pQ, 11, 20, 12 



9 The transform back from k to fco was n °t effected in the derivation presented in leading to a missing 
7|| -factor in the CP-violating force. The gamma-factor was included in JjiJ however. 
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5 Applications to supersymmetric models 



When studying electroweak baryogenesis in supersymmetric models an important source 
comes from CP-violating dynamics of charginos in presence of a bubble wall. In this case 
a formalism that describes the dynamics of mixing fermions is required. Introducing the 
mixing does not bring any extra complications w.r.t. to extension to 3+1 dimensions and 
hence, proceeding just as in KPSW-I we find the semiclassical Boltzmann equation 

d t fsi± + Vsi± ■ dxfsi± + F si ±d k J si ± = , (62) 

where f S i + and / S j_ are the distribution functions for fermions and antifermions of spin s, 
respectively, written in the propagating basis in which the mass matrix is diagonal. As it was 
the case in Eq. ([38|), / S j_'s are assumed to satisfy the condition f si ± = f S i±{k^ z,t — vu ■ x\\). 

The group velocity of the i'th quasiparticle eigenstate equals v s i± = k/u S i± and the 
CP-violating semiclassical force reads 



sii: 



2u si± 2^(^-^)1/2 



where |Mj| 2 B^ is a diagonal matrix which can be computed from the original mass matrix 
and the rotation matrix U diagonalizing MM^ as follows: 

|M d | 2 9' = \m{UMM l] U ] ) d . (64) 

For example, in the case of charginos the mass term reads 

^ R M ^ L + h.c. , (65) 

where = {W^,hi R ) T and = {Wi,h2 L ) T are chiral fields consisting of winos and 
higgsinos. The mass matrix is conveniently parametrized as 

where Hi = hie tdl and H 2 = h 2 e ie2 are the Higgs field vacuum expectation values and \x and 
m2 are the soft supersymmetry breaking parameters. For a realistic choice of parameters 
there is no spontaneous CP-violation in the MSSM, so to good approximation we can take 
the Higgs vev's to be real |22], ^3| . In this case [^(| f| U 

\ Ml \ 2 Q[ m 2 ± e' ± = T^(fxm 2 )(hih 2 y , (67) 
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where 9+ (0_) corresponds to the higgsino-like state when |/i| > |m 2 | (|/x| < |m 2 |), and 
A = m 2 + — vn?_ is the mass splitting of the chargino mass eig enstates @ g, g. CP violat ion 
is here mediated via the parameters /x, m 2 and may be large [f24}| - Yet another interesting 
CP-violating source arises when the Higgs condensate develops a CP-violating phase, as it 
is in the case of the NMSSM |12| (for details of computation see H). 



6 Conclusions and discussion 

We have presented a first principle derivation of the semiclassical limit of transport equations 
for a collisionless gas of fermions in a spatially varying background. Our analysis is based 
on the gradient expansion in the relative coordinate of the Kadanoff-Baym equations of 
motion |J for the out-of-equilibrium two point function G < . The most crucial step in our 
derivation was to show that to first order in H the Wigner transformed function G K has 
a spectral solution. The transport equations for the corresponding on-shell excitations are 
shown to reduce to the usual Liouville form, and contain a CP-violating semiclassical force 



of order K. We also showed how the same force can be obtained by WKB methods |L0 , [LT 



20| , H [12| . The force is suitable for baryogenesis calculations based on the charge transport 



picture |25| at a first order electroweak phase transition. 

The present work extends our earlier 1+1 dimensional results [|3| to the more realistic 
case of a planar phase transition wall in 3+1 dimensions, relevant for EWBG. The problem 
becomes analytically tractable only upon realising that the spin orthogonal to the moving 
planar wall remains a good quantum number for particles moving in arbitrary directions with 
respect to the wall. Thanks to this conserved quantity the problem can be reduced to the 
dynamics of one spin dependent distribution function for fermions, and one for antifermions. 
We assumed planar symmetry and stationarity appropriate for EWBG. For completeness 
we also discussed how our formalism can be extended to include particular time and space 
dependent transients which conserve spin. 

The spectral solution for G < can only be found up to order h. Beyond order h the 
spectral decomposition ansatz is not consistent and the transport equations do not have the 
semiclassical Liouville limit. It is therefore very fortunate that in at least some of the most 
interesting cases || the phase transition walls are rather wide £ wa u 3> ^dB, where is the 
de Broglie wave length ~ 1/T of a particle in the plasma, and so the semiclassical equation 
is expected to be a reasonable approximation. 
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The case of mixing fermions is relevant for example in supersymmetric baryogenesis sce- 
narios, in which baryogenesis can be mediated by mixing of charginos or neutralinos and 
other species. Supersymmetric models, in contrast to the minimal version of the standard 
model J26|, are shown to be viable candidates for EWBG fl2"7 |. We therefore derived ex- 



plicitly the CP-violating force for a general case of mixing fermions and for charginos in 
particular. With this, the present work covers the proof of the semiclassical methods used 
for baryogenesis calculations in the literature so far |9|, |12| . The directions for further work 
include the inclusion of the self-consistent gauge (hyperelectric) fields || in 3+1 dimensions 
and the study of interactions, which introduce new collisional baryogenesis sources 0. 
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